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Stochasticity effects on derivation of physical 
parameters of unresolved star clusters 
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Abstract. We developed a method for a fast modeling of broad-band UBVRI integrated 
magnitudes of unresolved star clusters and used it to derive their physical parameters (age, 
mass, and extinction). The method was applied on M33 galaxy cluster sample and con- 
sistency of ages and masses derived from unresolved observations with the values derived 
from resolved stellar photometry was demonstrated. We found that interstellar extinction 
causes minor age-extinction degeneracy for the studied sample due to a narrow extinction 
range in M33. 
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1. Introduction 

Star clusters are essential objects to constrain 
star formation history of their host galaxy. 
Traditional way to derive physical parameters 
of unresolved clusters is to compare their inte- 
grated broad-band photometry magnitudes to 
the magnitudes of simple stellar population 
models. However, the accuracy of parameter 
derivation is limited due to degeneracy be- 
tween the parameters and stochastic sampling 
of stars in clusters. 

Recently methods, which take stochasticity 
into account, for computing integrated cluster 
magnitudes by using Monte-Carlo sampling of 
stellar initial mass function, were introduced 
by Popescu & Hanson (2010) and Fouesneau 
& Lancon (2010). A large grid of cluster mod- 
els was built to cover wide ranges of age, 
mass, and extinction (metallicity assumed to be 
known). 



Physical parameters of clusters are then de- 
rived by comparing observations to the model 
grid. The first method finds a model which 
best matches observations, while the second 
is more accurate but much slower, because 
it builds probability maps in the age-mass- 
extinction space by exploring all the nodes of 
the grid and selects the most probable parame- 
ters for a given observation. 



We present a new method (de Meulenaer 
et al., 2013) that also builds probability maps 
of parameters, but does not require to explore 
all the nodes of the model grid, hence com- 
putation time is reduced significantly. We ap- 
plied this method on a cluster sample of M33 
galaxy studied by San Roman et al. (2009), Ma 
(2012), and new photometry from Ma (2013). 
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Fig. 1. Scheme of the method to derive physical pa- 
rameters of star cluster; see the text for description. 



2. Method to derive parameters 

The method is sketched in Fig. 1. It is fast, be- 
cause we do not explore all the models of the 
grid in photometric parameter space (panel a), 
but restrict to the ones located close to the ob- 
served absolute magnitudes of cluster, shown 
by circle, the size of which depends on photo- 
metric errors, cr. 

We use Padova isochrones by Girardi et al. 
(2010) to build cluster models. The model grid 
(panel b) has fixed metallicity (Z = 0.008), and 
contains: 71 steps in age, log(£/yr) = 6.6^-10.1; 
61 in mass, log(m/M©) = 2-^5; and 51 in ex- 
tinction, E(B - V) - -r 1. Each node of the 
grid (panel b) contains 1000 models. 

Panel (a) shows fraction of photometric pa- 
rameter space (without extinction) in plane of 
Mu vs Mb, with cluster's observation and error 
bars. All models located within the 3-cr area, 
which corresponds to nodes shown in panel 
(b), are selected. Distributions of mass and age 
shown in panels (c, d) are derived from the se- 
lected models. Most probable values are shown 
by solid vertical lines with confidence intervals 
(dashed lines). 

de Meulenaer et al. (2013) performed test 
on a sample of 10 4 artificial clusters with 
Gaussian photometric errors of 0.05 mag ap- 
plied for each passband to quantify accuracy of 
the method. Fig. 2 shows, as gray background, 



age (panels c, e), mass (panels d, f), and ex- 
tinction (panel b) derived for artificial sam- 
ple. The two streaks above and below the one- 
to-one relation (panel c) in the range of 8 < 
log(?/yr) < 9.5, first reported by Fouesneau & 
Lancon (2010), are a sign of the age-extinction 
degeneracy. 

3. Application on M33 clusters 

San Roman et al. (2009) derived age, mass, 
and extinction of 161 clusters in M33 galaxy 
using HST by isochrone fit to CMDs of their 
resolved stellar photometry. Ma (2012) pro- 
vided UBVRI broad-band integrated photom- 
etry for 40 of them. Recently, Ma (2013) ex- 
tended sample by 101 additional clusters and 
here we present results following study by de 
Meulenaer et al. (2013). Sample is presented in 
Fig. 2a; note that B band magnitudes of clusters 
were corrected by 0.1 mag. 

We adopted the M33 distance modulus of 
(m - Af)o = 24.54 (McConnachie et al, 2005). 
Extinction range of the M33 galaxy is assumed 
to be between the foreground Galactic line-of- 
sight extinction of E(B - V) = 0.04 (Schlegel 
et al., 1998) and 0.3, derived from supergiants 
(U et al., 2009) and HII regions (Rosolowsky 
& Simon, 2008). Interstellar extinction law of 
LMC (Gordon et al, 2003) is used to red- 
den the models. Note, that extinction range for 
clusters derived by San Roman et al. (2009) is 
more narrow (Fig. 2b) than expected in M33. 

The M33 cluster parameters were derived 
using two different allowed extinction ranges 
for the model grid: 1) narrow, E(B-V) = [0.04, 
0.30], and 2) wide, E(B - V) = [0.04, 1.0], as 
one used for the M31 case by de Meulenaer 
et al. (2013). Panels (b, c, d) of Fig. 2 show the 
extinction, age and mass derived for the narrow 
extinction range. Panels (e, f) show parameters 
derived with the wide range. Our results, com- 
pared to the values based on isochrone fit to 
CMDs by San Roman et al. (2009), show con- 
sistency. 

Differences are relatively small between 
the two kinds of parameter solutions using nar- 
row and wide extinction ranges of the model 
grid; panels: (c) vs (e), (d) vs (f). This sug- 
gests that the age-extinction degeneracy plays 
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Fig. 2. M33 galaxy star clusters common to catalogs of (San Roman et al., 2009, S09) and: 1) Ma (2012) 
(40 stars), 2) Ma (2013) (84 circles). Panel (a): model grid (without extinction) in the background; reddening 
vector follows LMC extinction law; line traces Padova SSP of Z = 0.008. Extinction (b), age (c, e), and 
mass (d, f) derived for clusters using integrated photometry of Ma (2012, 2013), compared to the values 
based on resolved stellar photometry by S09. Panels (b, c, d) are the case of narrow allowed extinction 
range, E(B - V) = [0.04,0.3]. Panels (e, f) - of wide range, E(B - V) = [0.04, 1.0]. Background shows 
results (derived vs true parameters) for artificial cluster sample; see text for description. 
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a minor role in this M33 cluster sample. 
Conversely, de Meulenaer et al. (2013) showed 
that degeneracy plays a major role in the M31 
galaxy preventing accurate parameter deriva- 
tion for a large fraction of clusters. 
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4. Conclusions 

The presented method takes into account 
stochastic sampling of stars in unresolved star 
clusters when deriving their physical parame- 
ters (age, mass, and extinction) using broad- 
band photometry. The M33 galaxy star cluster 
sample studied by San Roman et al. (2009) us- 
ing resolved stellar photometry (isochrone fit 
to CMDs), was analyzed based on the UBVRI 
integrated photometry from Ma (2012, 2013), 
and shows consistent ages and masses. 

By studying the cluster sample succes- 
sively in narrow and wide allowed extinction 
ranges, we demonstrate that derived cluster pa- 
rameters indicate low sensitivity to the age- 
extinction degeneracy, contrarily to what is 
the case for the M31 galaxy star clusters (de 
Meulenaer et al., 2013). 
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